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Abstract 

We discuss the statistical mechanics of a system of self-gravitating fermions in a space 
of dimension D. We plot the caloric curves of the self-gravitating Fermi gas giving the 
temperature as a function of energy and investigate the nature of phase transitions as a 
function of the dimension of space. We consider stable states (global entropy maxima) as 
well as metastable states (local entropy maxima). We show that for D > 4, there exists 
a critical temperature (for sufficiently large systems) and a critical energy below which 
the system cannot be found in statistical equilibrium. Therefore, for D > 4, quantum 
mechanics cannot stabilize matter against gravitational collapse. This is similar to a 
result found by Ehrenfest (1917) at the atomic level for Coulombian forces. This makes 
the dimension D = 3 of our universe very particular with possible implications regarding 
the anthropic principle. Our study enters in a long tradition of scientific and philosophical 
papers who studied how the dimension of space affects the laws of physics. 

1 Introduction 

The statistical mechanics of systems with long-range interactions is currently a topic of ac- 
tive research [Ij. Among long-range interactions, the gravitational force plays a fundamental 
role. Therefore, the developement of a statistical mechanics for self-gravitating systems is of 
considerable interest [2]. In this context, a system of self-gravitating fermions enclosed within 
a box provides an interesting model which can be studied in great detail [31 H]. This model 
incorporates an effective small-scale cut-off played by the Pauli exclusion principle and a large 
scale cut-off played by the confining box (other forms of confinement could also be considered). 
The statistical mechanics of this system is rigorously justified and presents a lot of interesting 
features which are of interest in statistical mechanics [5] and astrophysics Its detailed study 
is therefore important at a conceptual and practical level. 

In a preceding paper [1], we have discussed the nature of phase transitions in the self- 
gravitating Fermi gas in a space of dimension D = 3. Our study was performed in both 
microcanonical and canonical ensembles and considered an arbitrary degree of degeneracy rel- 
ative to the system size. This study completes previous investigations by Hertel & Thirring 
[3] who worked in the canonical ensemble and considered small system sizes. At high temper- 
atures and high energies, the system is in a gaseous phase and quantum effects are completely 
negligible. At some transition temperature or transition energy (for sufficiently large 
system sizes), a first order phase transition is expected to occur and drive the system towards 
a condensed phase. However, gaseous states are still metastable below this transition point 
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and gravitational collapse will rather occur at a smaller critical temperature Tc (Jeans tem- 
perature) [7] or critical energy Ec (Antonov energy) [HI [9l [2] at which the metastable branch 
disappears (spinodal point). The end-state of the collapse is a compact object with a "core- 
halo" structure. Typically, it consists of a degenerate nucleus surrounded by a "vapour" . The 
nucleus (condensate) resembles a white dwarf star [TU]. At non-zero temperature, this compact 
object is surrounded by a dilute atmosphere. Therefore, when quantum mechanics is properly 
accounted for, there exists an equilibrium state (global maximum of entropy or free energy) 
for each value of accessible energy and temperature. The condensate results from the balance 
between gravitational contraction and quantum pressure. As first noticed by Fowler in his 
classical theory of white dwarf stars, quantum mechanics is able to stabilize matter against 
gravitational collapse. 

One object of this paper is to show that this conclusion is no more valid in a space of 
dimension D > A. For a system of mass M enclosed within a box of radius R, there exists a 
critical temperature (for sufficiently large R) and a critical energy below which the system can- 
not be found at statistical equilibrium. This is like the Antonov instability for self-gravitating 
classical particles in D = 3 [HI El 12] but it now occurs for fermions. Therefore, quantum me- 
chanics cannot arrest gravitational collapse in D > A. This result is connected to our previous 
observation [12j that a classical white dwarf star (a polytrope of index 713/2 = D/2) becomes 
unstable for D > 4. Interestingly, this result is similar to that of Ehrenfest [13] who consid- 
ered the stability of atomic structures (in Bohr's model) for different dimensions of space and 
concludes that D < A is required for stability. In this paper, we determine the caloric curve of 
the self-gravitating Fermi gas for an arbitrary dimension of space and an arbitrary degree of 
degeneracy (or system size). We exhibit particular dimensions that play a special role in the 
problem. The dimension D = 2 is critical because the results established for D ^ 2 cannot be 
directly extended to D = 2 [T3|. Furthermore, in D = 2 the radius of a white dwarf star is in- 
dependent on its mass and given in terms of fundamental constants hy R = 0.27 hm~^^'^G^^^^. 
The dimension D = 4 is also critical because it is the dimension at which classical white dwarf 
stars become unstable. At this particular dimension, their mass is independent on radius and 
can be expressed in terms of fundamental constants as M = 1.44 10~^/i^m~^G~^. Mathemati- 
cally, this is similar to Chandrasekhar's limiting mass [15] for relativistic white dwarf stars in 
D = 3. The dimension D = 2(1 + y^) is also particular because at this dimension, the white 
dwarf stars cease to be self-confined and have infinite mass. Finally, Z) = 10 is the dimension at 
which the caloric curve of classical isothermal spheres loses its characteristic spiral nature [T^ . 
Although we systematically explore all dimensions of space in order to have a complete picture 
of the problem, only dimensions D = 1, D = 2 and D = 3 are a priori of physical interest. 
The dimension D = 1 is considered in cosmology and in connexion with shell models, and the 
dimension D = 2 can be useful to describe filaments or ring structures with high aspect ratio. 
Two-dimensional gravity is also of interest for its properties of conformal invariance and for its 
relation with two-dimensional turbulence [5|. Non-integer dimensions can arise if the system 
has a fractal nature. 

The paper is organized as follows. In Sec. [2], we determine the thermodynamical parameters 
of the self-gravitating Fermi gas in dimension D. The Fermi-Dirac entropy is introduced from a 
combinatorial analysis. In Sec. [3l we consider asymptotic limits corresponding to the classical 
self-gravitating gas and to completely degenerate structures (white dwarfs). We emphasize the 
importance of metastable states in astrophysics and explain how they can be taken into account 
in the theory (see also [16]). We also discuss the thermodynamic limit of the self-gravitating 
quantum gas and compare it with the thermodynamic limit of the self-gravitating classical gas 
in the dilute limit [T7]. In Sec. HI we provide a gallery of caloric curves of the self-gravitating 
Fermi gas in different dimensions of space. Rigorous mathematical results on the existence of 
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solutions of the Fermi- Poisson equation have been obtained by Stahczy [IB]. Finally, in the 
conclusion, we place our study in a more general perspective. We give a short historical account 
of scientific and philosophical papers who studied the role played by the dimension of space in 
determining the form of the laws of physics. These works tend to indicate that the dimension 
D = 3 oi our universe is very particular. This is also the result that we reach in our study. 
These remarks can have implications regarding the anthropic principle. 



2 Thermodynamics of self-gravitating i^-fermions 
2.1 The Fermi-Dirac distribution 

We consider a system of N fermions interacting via Newtonian gravity in a space of dimension 
D. We assume that the mass of the configuration is sufficiently small so as to ignore relativistic 
effects. Let /(r, v,t) denote the distribution function of the system, i.e. f{r,\',t)d^rd^\ gives 
the mass of particles whose position and velocity are in the cell (r, v; r + d^r, v + d^v) at time 
t. The integral of / over the velocity determines the spatial density 

p = J fd^v, (1) 

and the total mass of the configuration is given by 

M = [ pd^r, (2) 



where the integral extends over the entire domain. On the other hand, in the meanfield ap- 
proximation, the total energy of the system can be expressed as 

E = ^J fv^d^rd^v + i y pM^r = K + W, (3) 

where K is the kinetic energy and W the potential energy. The gravitational potential $ is 
related to the density by the Newton-Poisson equation 

A$ = SoGp, (4) 

where Sd = 2tx^/'^ /V{D /2) is the surface of a unit sphere in a space of dimension D and G is 
the constant of gravity (which depends on the dimension of space). 

We now wish to determine the most probable distribution of self-gravitating fermions at 
statistical equilibrium. To that purpose, we divide the individual phase space {r, v} into a very 
large number of microcells with size {h/m)^ where h is the Planck constant (the mass m of 
the particles arises because we use v instead of p as a phase space coordinate). A microcell 
is occupied either by or 1 fermion (or g = 2s + 1 fermions if we account for the spin). We 
shall now group these microcells into macrocells each of which contains many microcells but 
remains nevertheless small compared to the phase-space extension of the whole system. We 
call u the number of microcells in a macrocell. Consider the configuration {rij} where there 
are ni fermions in the 1*^* macrocell, n2 in the 2'^'^ macrocell etc., each occupying one of the 
i> microcells with no cohabitation. The number of ways of assigning a microcell to the first 
element of a macrocell is u, to the second z/ — 1 etc. Since the particles are indistinguishable, 
the number of ways of assigning microcells to all particles in a macrocell is thus 

1 u\ 
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To obtain the number of microstates corresponding to the macrostate {rij} defined by the 
number of fermions Ui in each macrocell (irrespective of their precise position in the cell), 
we need to take the product of terms such as ([5]) over all macrocells. Thus, the number of 
microstates corresponding to the macrostate {nj}, i.e. the probability of the state {nj}, is 



W{{n,}) = X{- 



(6) 



ni\{u -ni)\' 

This is the Fermi-Dirac statistics. As is customary, we define the entropy of the state {rii} by 

S{{n,})=\nW{{n,}). (7) 

It is convenient here to return to a representation in terms of the distribution function giving 
the phase-space density in the i-th macrocell 



fi = f{r^, Vj 



rii m 



V 



where we have defined r/o = vrP^^ jh^ ^ which represents the maximum value of / due to Pauli's 
exclusion principle. Now, using the Stirling formula, we have 



lnVr(K}) ~ ^z/(ln, 
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Passing to the continuum limit u 
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0, we obtain the usual expression of the Fermi-Dirac 



1 - 



In 1 



^0 V VoJ V VoJ ) (^)^ 

If we take into account the spin of the particles, the above expression remains valid but the 
maximum value of the distribution function is now tjq = gm^~^^/h^, where = 2s -|- 1 is the 
spin multiplicity of the quantum states (the phase space element has also to be multiplied by 
g). An expression of entropy similar to (fTOl) . but arising for a completely different reason, 
has been introduced by Lynden-Bell in the context of the violent relaxation of collisionless 
stellar systems [191 1201 121]. In that context, tiq represents the maximum value of the initial 
distribution function and the actual distribution function (coarse-grained) must always satisfy 
f Vo virtue of the Liouville theorem. This is the origin of the "effective" exclusion 
principle in Lynden-Bell's theory, which has nothing to do with quantum mechanics. Since the 
particles (stars) are distinguishable classical objects (but subject to an exclusion principle in 
the collisionless regime), Lynden-Bell's statistics corresponds to a 4-th form of statistics (in 
addition to the Maxwell-Boltzmann, Fermi-Dirac and Bose- Einstein statistics). However, for a 
single type of phase element rjo, it leads to the same results as the Fermi-Dirac statistics. We 
also recall that in the non-degenerate (or classical) limit / -C r/o, the Fermi-Dirac entropy (fTO!) 
reduces to the Boltzmann entropy 
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Now that the entropy has been precisely justified, the statistical equilibrium state (most 
probable state) of self-gravitating fermions is obtained by maximizing the Fermi-Dirac entropy 
( JTOl) at fixed mass ([2]) and energy ([3]): 



Max S[f] I E[f]=E,M[f]=M. 



(12) 
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Introducing Lagrange multipliers 1/T (inverse temperature) and /i (chemical potential) to sat- 
isfy these constraints, and writing the variational principle in the form 

6S-^6E + ^6N = 0, (13) 

we find that the critical points of entropy correspond to the Fermi-Dirac distribution 

/ = (14) 

where A = e~^'^ is a strictly positive constant (inverse fugacity) and /? = is the inverse 
temperature. Clearly, the distribution function satisfies / < r/o, which is a consequence of 
Pauli's exclusion principle. 

So far, we have assumed that the system is isolated so that the energy is conserved. If now 
the system is in contact with a thermal bath (e.g., a radiation background) fixing the temper- 
ature, the statistical equilibrium state minimizes the free energy F = E — TS, or maximizes 
the Massieu function J = S — PE, at fixed mass and temperature: 

Max J[f] I M[f] = M. (15) 

Introducing Lagrange multipliers and writing the variational principle in the form 

6J + ^6N = 0, (16) 

we find that the critical points of free energy are again given by the Fermi-Dirac distribution 
( IT4ll . Therefore, the critical points (first variations) of the variational problems ( fT2l) and (fTSl) 
are the same. However, the stability of the system (regarding the second variations) can 
be different in microcanonical and canonical ensembles. When this happens, we speak of a 
situation of ensemble inequivalence The stability of the system can be determined by a 
graphical construction, by simply plotting the caloric curve/series of equilibria l3{E) and using 
the turning point method of Katz |22i i23j . 



2.2 Thermodynamical parameters 

Integrating the distribution function (fT^ over velocity, we find that the density of particles is 
related to the gravitational potential by 



„ c 9D/2-I 

ipm)^/^ -^i^/2-i(Ae^ ), (17) 

where denotes the Fermi integral 



We recall the identity 

m = ~In-i{t), {n>0), (19) 



which can be established from ( flSl) by an integration by parts. The gravitational potential is 
now obtained by substituting Eq. ( |T71) in the Poisson equation (jl]). We introduce the rescaled 
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distance ^ = [Sl2^''^-^G7]Q/{(3m)^l^-^Y/'^r and the variables ?/> = (3m{^-%) and k = Ae^™*o, 
where $o is the central potential. Thus, we get the D-dimensional Fermi-Poisson equation 



=^-/-(^^'"^)' (20) 



m = ^'(0) = 0. (21) 

As is well-known, self-gravitating systems at non-zero temperature have the tendency to 
evaporate. Therefore, there is no equilibrium state in a strict sense and the statistical me- 
chanics of self-gravitating systems is essentially an out-of-equilibrium problem. However, the 
evaporation rate is small in general and the system can be found in a quasi-equilibrium state 
for a relatively long time. In order to describe the thermodynamics of the self-gravitating Fermi 
gas rigorously, we shall use an artifice and enclose the system within a spherical box of radius 
R (the box typically represents the size of the cluster under consideration). In that case, the 
solution of Eq. (120!) is terminated by the box at the normalized radius 



a 



(22) 



(/3m)^/2-i 

For a spherically symmetric configuration, the Gauss theorem can be written 

_ GM{r) 
dr yD-i ' 



(23) 



where M{r) = pSof^ ^dr is the mass within the sphere of radius r. Applying this result at 
r = R and using the variables introduced previously we get 

^ = ^^^D^ = (24) 

This equation relates the dimensionless box radius a and the uniformizing variable k to the 
dimensionless inverse temperature r]. According to Eqs. (l22l) and (12^ . a and k are related to 
each other by the relation a'^r]^^'^^^ = fi or, explicitly. 



-D+2 , 2 

ao-2 7/';(a) =/iO-2, (25) 



p = 7]QJS%2D-^GDMD-m'^i*-^), (26) 



where 



is the degeneracy parameter [20]. It should not be confused with the chemical potential. We 
shall give a physical interpretation of this parameter in Sec. 14. 2[ 

The calculation of the energy is a little more involved. First, we introduce the local pressure 

p = 1 / /t;2rf^v. (27) 



D „ 

Using the Fermi-Dirac distribution function f|T^ . we find that 



The kinetic energy K = (-D/2) J pd^r can thus be written 



P = ^?P^lDn{ke^)- (28) 



KR^-^ a D-2 



GM^ 



r lD/2{ke'i'-^^^)e-'di. (29) 
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In order to determine the potential energy, we use the D-dimensional version of the Virial 
theorem [12]. For D 7^ 2, it reads 



2K + (D -2)W = DVoR^piR), (30) 

where Vd = Sd/D is the volume of a hypersphere with unit radius (the case D = 2 will be 
considered specifically in Sec. 14. 6p . Using the expression of the pressure f l28!) at the box radius 
R, we get 

WR''-' _ 2 a— ^(,) _ 2KR^-' 

GM^ ~ D{D-2) > (Z}-2)GM2- ^^^> 

Combining Eqs. ( !29l) and ( !3TI) . we finally obtain 

(32) 

For D = 3, Eqs. (|2^ and fl32|) return the expressions derived in [201 II]- For a given value 
of fi and k, we can solve the ordinary differential equation fl20|) until the value of a at which 
the condition fl25]) is satisfied. Then, Eqs. (1211) and fl52]) determine the temperature and the 
energy of the configuration. By varying the parameter k (for a fixed value of the degeneracy 
parameter /x), we can determine the full caloric curve/series of equilibria j3{E). Extending the 
results of [1] in D dimensions, the entropy of each configuration, parameterized by a, is given 
by 

2D 

jk - -'-W^^^ ^ ^.(") ^ ^ ^ ^ - |^^/..(-«"'). (33) 
and the free energy by 

F = E-TS. (34) 

In the microcanonical ensemble, a solution is stable if it corresponds to a maximum of en- 
tropy S[f] at fixed mass and energy. In the canonical ensemble, the condition of stability 
requires that the solution be a minimum of free energy F[f] at fixed mass and temperature. 
This meanfield approach is exact in a thermodynamical limit such that +00 with /i, 

r/, A fixed. If we fix r/o (i.e. h) and G, this implies that RN(D-2)/{D{i-D)) ^ y^_4/(D(4-D))^ 
^jY-(4D-d2+4)/(d(4-d))^ jjY~i approach a constant value for ^ +00 (the free 

energy F scales as A^(*^-^~^^+^)/(-^(^~^))). This is the quantum thermodynamic limit (QTL) 
for the self-gravitating gas [U |23]. The usual thermodynamic limit N,R —>■ +00 with N/R^ 
constant is clearly not relevant for inhomogeneous systems whose energy is non-additive. 



3 Asymptotic limits 

3.1 The non degenerate limit (/i = 00) 

Before considering the case of an arbitrary degree of degeneracy, it may be useful to discuss 
first the non degenerate limit corresponding to a classical isothermal gas (^ — > 0). For / -C r^o, 
the distribution function (IT^ reduces to the Maxwell-Boltzmann formula 

/ = ^e-^'"(^+*), (35) 
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which can be written more conveniently as 



/=(^)"V)e-/^-4. (36) 
The density profile can be written 

P = Poe-^(«\ (37) 
where po is the central density, ^ is the normalized distance 

^ = {SoGPrnpoY/'r, (38) 

and ip is the solution of the D-dimensional Emden equation 

1 d f^D-^d_^\^^^,^ (39) 



with boundary conditions 

^(0) = ip\0) = 0. (40) 

This equation can be obtained from Eq. ( !20l) by taking the limit k +oo and using the 
limiting form of the Fermi integral 

/„(t) ~ ^r(n + i), (t-^+oo). (41) 

From Eq. (136!) . we check that the local equation of state of a classical self-gravitating isothermal 
gas is p{r) = ^^ksT whatever the dimension of space. The thermodynamical parameters are 
given by 

T] = a^'{a), (42) 

D{A~D) 1 1 e-'^^") 

2{D - 2) aip'{a) ~ D-2ip'{a) 

° ~ lnr/-21na + ^(a) + — 2Ar/, (44) 



Nks 2 ' ' D-2 

^0 , , /27r^/2\ ^ D 



Nk. ^''^ 

where a = {SoGPrnpoY^'^R is the normalized box radius. For D = 2, the thermodynamical 
parameters can be calculated analytically [H]. Introducing the pressure at the box P = p{R), 
the global equation of state of the self-gravitating gas can be written 

^ = l^e-*(">. (46) 



NkeT D 7] 

We recall that the foregoing expressions can be expressed in terms of the value of the Milne 
variables Uq = u{a) and Vq = v{a) at the normalized box radius [3 [25] • The structure and 
the stability of classical isothermal spheres in D dimensions have been studied in detail in [13] . 
The classical thermodynamic limit (CTL) of self-gravitating systems, or dilute limit [17], is 
such that +oo with i], A fixed. If we take /? ~ 1, this implies that R ~ N'^/i^-'^) and 

E, S, J,F ~ A^. The physical distinction between the quantum thermodynamic limit (QTL) 
and the classical thermodynamic limit (CTL) is related to the existence of long-lived gaseous 
metastable states as discussed in 
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3.2 The completely degenerate limit 

For P — > +00 (i.e., T = 0), the distribution function (fT^ reduces to a step function: / = 779 if 
V < vp and / = if f > fi?, where f_F(r) = yjH^pjm — $) is the local Fermi velocity. In that 
case, the density and the pressure can be explicitly evaluated: 

^^Sj.v^'-Hv = VoSd^, (47) 







1 /■''^ S'r. 7;-°+2 



"'dL ""S^^-V^ = „„^^. (48) 

Eliminating the Fermi velocity between these two expressions, we find that the equation of 
state of a cold Fermi gas in D dimensions is 

This equation of state describes a £)-dimensional classical white dwarf star (throughout this pa- 
per, we shall call "white dwarf star", or "fermion ball", a completely degenerate self-gravitating 
system. This terminology will be extended to any dimension of space). In D = 3, classical 
white dwarf stars are equivalent to polytropes with index n = 3/2 [11]. In D dimensions, 
classical "white dwarf stars" are equivalent to polytropes with index 



D , . 

^3/2 = —■ (50) 

The structure and the stability of polytropic spheres in D dimensions have been studied in detail 
in p^. It is shown that a polytrope of index n is self-confined for n < = {D + 2)/{D — 2) and 
stable for n < = D/{D — 2). Therefore, white dwarf stars (n = 11^/2 = D/2) are self-confined 
only for D < 2(1 + -\/2) and they are stable only for D < A. For D > A, quantum mechanics is 
not able to stabilize matter against gravitational collapse. Thus, D = 4 is a critical dimension 
regarding gravitational collapse. D = 2 is also critical [11]. Therefore, the dimension of space 
of our universe 2 < D = 3 < 4 lies between two critical dimensions. 

We now introduce dimensionless parameters associated with ^3/2 polytropes which will be 
useful in the sequel. Their density profile can be written 

p(r)=po^^/'(0, (51) 
where po is the central density, ^ is the normalized distance 

2SdGpq ° 



K{D + 2) 

and 6 is solution of the Z)-dimensional Lane-Emden equation 



r, (52) 



1 d f^^^,de\_^^^,^ ^^^^ 



with boundary conditions 



^(0) = 1, e'{0) = 0. (54) 
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This equation can be obtained from Eq. fl2UI) by taking the hmit k —>■ and using the hmiting 
form of the Fermi integral 

(-lnt)"+i , , 

~ HtTi ' ^ ^^^^ 

For D < 2(1 + \/2), the solution of the Lane-Emden equation (!53l) vanishes at a finite distance 
^1 defining the radius R^, of the white dwarf star {complete polytrope). Using the results of [12], 
the mass-radius relation of D-dimensional white dwarf stars is given by 

where we have defined 

^D/2 = -eP^'te)- (57) 

For 2 < D < 4, the mass M decreases with the radius i?* while for D < 2 and for A < D < 
2(1 + V2) it increases with the radius (see Fig. [1]). The mass-radius relation ( l56l) exhibits the 
two critical dimensions of space D = 2 and D = A discussed previously. For D = 2, the radius 
is independent on mass and for D = 4, the mass is independent on radius (see Sec. H]). The 
energy of a self-confined white dwarf star is 

E = -Ad/2— dZ2, (58) 
-fx* 

where 

^^/^ - (D-2)(4 + 4D-D2)- ^^^^ 

We note that the energy of a white dwarf star vanishes for D = A. According to Poincare's 
theorem [10], this determines the onset of instability. We thus recover the fact that complete 
white dwarf stars are unstable for D > 4 [12|. 

For D > 2(1 + V2), the density of a 123/2 polytrope never vanishes (as ^3/2 > ^5) and we 
need to confine the system within a box of radius R {incomplete polytrope) to avoid the infinite 
mass problem. In that case, the white dwarf star exerts a pressure against the box. White 
dwarf stars with R^, > R when D < 2(1 -|- \/2) are also incomplete. They are arrested by the 
box at the normalized radius = a. with a = {2SdGp~^ '^'^^^ /[K{D + 2)]}^/^i?. As shown in 
[12] . the normalized mass and the normalized energy of the configuration parameterized by a 
are given by 



M 

Vp = ^ 

Do 



2S £)G 



D 

D-5 



1 



D+2 



A = 



[K{D + 2)\ 
ER^-^ -2 (D{4-D 



-aT^e'{a), (60) 



GM^ D'^ -AD -4:{2{D -2 



l + (fl-2)-''<") 



a9'(a) 



2~De{a 
2 + D e'{a) 



D+2 
2 



(61) 



In the present context, the normalized mass rjp is related to the degeneracy parameter /i by 
the relation ^ 

2/i\ ^ 



On the other hand, using Eqs. (jSS]) and (1581) . the normalized mass and the normalized energy 
of a self-confined white dwarf star with R^ < R (complete polytrope) are given by 

(D-4)D 

RA^^ 



VP = { ) (63) 
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R/R„ 



Figure 1: The mass-radius relation for complete white dwarf stars (T = 0) in different di- 
mensions of space. It clearly shows that the dimension D = 3 in surrounded by two critical 
dimensions D = 2 and D = 4 at which either the radius or the mass is constant. 



A = Xd/2 i^j^J . (64) 
Eliminating i?^, between these two relations, we obtain the "mass-energy" relation 



Av^"' =>^D/2iu;D/2)'^, (65) 
which will be useful in our subsequent analysis. 



4 Caloric curves in various dimensions 
4.1 Series of equilibria and metastable states 

We shall now determine the caloric curve (3{E) of the self-gravitating Fermi gas as a function 
of the degeneracy parameter /i for any dimension of space D. This study has already been 
performed for Z) = 3 in i4j. The critical points of the Fermi-Dirac entropy S[f] at fixed E and 
M (i.e., the distribution functions /(r, v) which cancel the first order variations of 5* at fixed E, 
M) form a series of equilibria parameterized by the uniformizing variable k. At each point in 
the series of equilibria corresponds a temperature j3 and an energy E determined by Eqs. flMl) 
and fl32|) . In this approach, (3 is the Lagrange multiplier associated with the conservation of 
energy in the variational problem f|T3|) . It has also the interpretation of a kinetic temperature 
in the Fermi-Dirac distribution f[T^ . We can thus plot P{E) along the series of equilibria. 
There can be several values of temperature (3 for the same energy E because the variational 
problem (fT2l) can have several solutions: a local entropy maximum (metastable state), a global 
entropy maximum, and one or several saddle points. We shall represent all these solutions on 
the caloric curve because local entropy maxima (metastable states) are in general more physical 
than global entropy maxima for the timescales achieved in astrophysics. Indeed, the system 
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can remain frozen in a metastable gaseous phase for a very long time. This is the case, in 
particular, for globular clusters and for the gaseous phase of fermionic matter (at high energy 
and high temperature). The time required for a metastable gaseous system to collapse is in 
general tremendously long and increases exponentially with the number of particles (thus, 
tiife +00 in the thermodynamic limit +00) [16j. This is due to the long-range nature 

of the gravitational potential. Therefore, metastable states are in reality stable states. At 
high temperatures and high energies, the global entropy maximum is not physically relevant 
[261 [271 [25l [21] . Condensed objects (e.g., planets, stars, white dwarfs, fermion balls,...) only 
form below a critical energy Ec (Antonov energy) |8l |9l [2] or below a critical temperature Tc 
(Jeans temperature) [7], when the gaseous metastable phase ceases to exist (spinodal point). 

4.2 The case 2 < L> < 4 

We start to describe the structure of the caloric curve of the self-gravitating Fermi gas for 
2 < D < A (specifically D = 3). Let us first consider the Fermi gas at T = (white dwarf 
stars). The A — rjp curve defined by Eqs. (pi . (IMl) and is represented in Fig. [21 In the 
present context, it gives the energy of the star as a function of its mass. Since the curve does 
not present turning points, all the white dwarf star configurations are stable. According to Eq. 
( 156|) . for 2 < D < 4, the mass M of a complete white dwarf star is a decreasing function of its 
radius i?^,. Therefore, if the system is enclosed within a box, there exists a characteristic mass 

M,{R) = R-^^ (66) 

such that for M > M^{R) the star is self-confined [R^ < R) and for M < M^[R), it is restricted 
by the box. In terms of the dimensionless mass r^p, complete polytropes correspond to 
'Hp ^ and incomplete polytropes to rjp < ujd/2- For 2 < D < 4, there exists a stable 
equilibrium at T = for all mass M. 

We now briefiy describe the caloric curve for arbitrary temperature and energy. A more 
complete description is given in |1] for D = 3. First, we note that, according to Eqs. (l26l) . ( l49i) 
and dMl), 

D(4.-D) 

where 

/i,(D) = ^(^B/2)'^. (68) 

Therefore, the degeneracy parameter /i can be seen as the ratio (with some power) between 
the size of the system R and the size i?* of a white dwarf star with mass M. Accordingly, a 
small value of corresponds to a large "effective" cut-off (played by Pauli's exclusion principle) 
or, equivalently, to a small system size. Alternatively, a large value of /i corresponds to a 
small "effective" cut-off or a large system size. This gives a physical interpretation to the 
degeneracy parameter. For yU +00 (i.e. /I ^ 0), we recover classical isothermal spheres. 
In that case, the caloric curve (3{E) forms a spiral. For finite values of /i, the spiral unwinds 
due to the infiuence of degeneracy and gives rise to a rich variety of caloric curves (Fig. [S]). 
For large systems, the caloric curve has a Z-shape ("dinosaur's neck") and for small systems 
it has a A^-shape. The phase transitions in the self-gravitating Fermi gas for D = 3 and the 
notion of metastable states, spinodal points, critical points, collapse, explosion, and hysteresis 
are discussed in [H [271 EH HE] . Similar notions are discussed in [28] for a hard spheres gas. 
The ground state of the self-gravitating Fermi gas (T = 0) corresponds to a white dwarf star 
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Figure 2: The mass-energy relation for white dwarf stars (T = 0) in D = 3. There exists 
an equihbrium state for all mass. The white dwarf star is self-confined if M > M^(i?) and 
box-confined if M < M^{R). 
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Figure 3: Caloric curve in D = 3 for different values of the degeneracy parameter (various 
system sizes). 
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Figure 4: Caloric curve in D = 3 for small values of the degeneracy parameter (small system 
sizes). For D < 4, there exists an equilibrium state for all temperatures T and all accessible 
energies E > Eground- 

configuration. For given fi, its structure (radius, energy) is determined by the intersection 
between the A — rjp curve in Fig. [2] and the line defined by Eq. fl62l) . The "white dwarf" is 
complete {R^: < R) for ^ > fi^{D) and incomplete (i?^, > R) otherwise. For > fJ.^:{D), the 
normalized energy of the white dwarf is given by 



This is the ground state of the self-gravitating Fermi gas corresponding to the asymptote in 
Fig. m (this asymptote exists for all curves in Fig. [3] but is outside the frame). For classical 
particles {h = 0), there is no equilibrium state if energy and temperature are below a critical 
threshold [Hl[9]- In that case, the system undergoes gravitational collapse and forms binaries 
(in microcanonical ensemble) or a Dirac peak (in canonical ensemble); see Appendices A and B 
of [H] and [25l [2^ [16] . For self-gravitating fermions, an equilibrium state exists for all values 
of temperature and for all accessible energies {E > Eground)- Gravitational collapse is arrested 
by quantum pressure as first realized by Fowler [TT]. We shall now show that this claim ceases 
to be true in dimension D > 4. 



We now consider the case 4 < D < 2(1 + \/2) (specifically D = 4.1). Let us first describe 
the Fermi gas at T = 0. The A — rjp curve defined by Eqs. (!60|) . (!6T|) and (165|1 is represented 
in Fig. [51 For D > 4, the curves ?7p(a) and A{a) associated to 113/2 polytropes have their 
extrema at the same point (see Appendix C of [12j). Therefore, the A-r]p curve presents a 
cusp at (Aq, r^p.c)- Past this point in the series of equilibria, 7x3/2 polytropes are unstable. 
According to Eq. fl56p . for D > 4, the radius i?^, of a self-confined white dwarf star increases 
with its mass. For M < (R) there exists self-confined white dwarf star configurations. In 
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(69) 



4.3 The case 4 < D < 2(1 + v^) 
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terms of the dimensionless mass rjp, this corresponds to rjp < ujd/2 (see Fig. Ej). However, 
such configurations are unstable since they he after the turning point [12]. Therefore, only 
incomplete (box confined) white dwarf stars can be stable in D > 4. Inspecting Fig. [5] again, 
we observe that these configurations exist only below a critical mass 



K{D + 2) 



D 
D-2 



(70) 



2S£)G 

For M > Mc{R), there is no equilibrium state at T = for D > 4. In terms of the dimensionless 
mass Tjp, equilibrium states exist only for rjp < rjp^dD)- 
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Figure 5: The mass-energy relation for white dwarf stars (T = 0) in 4 < D < 2(1 + v^) 
(specifically D = 4.1). Self-confined white dwarf stars are always unstable. Box-confined white 
dwarf stars exist only for M < Mc{R). For M > Mc{R), there is no equilibrium state. 

The caloric curve for arbitrary value of temperature and energy is represented in Fig. El 
For n — > +00, we recover the classical spiral [11]. For finite values of fi, there exists equilibrium 
solutions at all temperatures only if rjp < rjpc{D). Using Eq. (!62|) . this corresponds to 

/i < ^np,,{D)^ = fi.iD). (71) 



-A^GM^R^- 
^) below which 



If yU > jj,(.{D), or equivalently if M > M^{R), there exists a minimum energy = 
(which appears to be positive) and a minimum temperature = GM/ {r]cR^ 
there is no equilibrium state (the values of rjc and Ac depend on D and fi). In that case, 
the system is expected to collapse. This is similar to the Antonov instability (gravothermal 
catastrophe) for classical particles [3 [9]. Since we deal here with self-gravitating fermions, we 
could expect that quantum pressure would arrest the collapse. Our study shows that this is not 
the case for D > A. Quantum mechanics cannot stabilize matter against gravitational collapse 
anymore. 



4.4 The case L> = 4 

The dimension D = A is special because it is the dimension of space above which quantum 
pressure cannot balance gravity anymore. Therefore, D = 4 is critical and it deserves a partic- 
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Figure 6: Caloric curve in D = 4.1 for different values of the degeneracy parameter. For 
fi > fic{D), there is no equilibrium state if the temperature and the energy are too low. The 
reason for the "gap" at kc is explained in Fig. [71 




Figure 7: Graphical construction determining the value of a for given /i and k (in D = 4.1). 
According to Eq. ( l25i) . the normalized box radius a is solution of /(a) = 0, where /(^) = 

D+2 2 

C°~^Vk\0 ~ /^^"^- We see that a undergoes a discontinuity as A; kc. This gives rise to the 
"gap" in Fig. [6] for = 23. However, this gap is essentially a mathematical curiosity since the 
lower part of the curve (small k) is unstable anyway. 
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ular attention. First, consider the Fermi gas at T = 0. It corresponds to a polytrope of index 
^3/2 = '^3 [I2]. The A — rjp curve defined by Eqs. (!60|) . (!6T|) and (!65|) is represented in Fig. |8l 
Since the curve is monotonic the box-confined 71^/2 polytropes are stable and the complete 71,3/2 
polytropes are marginally stable. For D = 4, the mass of a self-confined white dwarf star is 
independent on its radius, see Eq. (13^ . It can be expressed in terms of fundamental constants 
as 

funru = -4^^^ ^ 1.44 10-2 (72) 



gSl m^G^ 



m 



where LJ2 — 11.2 (we have taken g = 2 in the numerical application). Mathematically, this is 
similar to Chandrasekhar's limiting mass for relativistic white dwarf stars equivalent to n = 3 
polytropes in D = 3 [15J. However, it is here a purely classical (i.e. nonrelativistic) result. 
Relativistic effects will be considered in a forthcoming paper [2S|. The energy of the self- 
confined white dwarf stars is £^ = 0. Considering Fig. [8] again, we see that incomplete white 
dwarf stars exist only for M < Mumu- In terms of the dimensionless mass 7]p, this corresponds 
to ?7 < rjpc = UJ2 — 11.2. For M > Mumu, there is no equilibrium state at T = 0. The 
caloric curve for arbitrary value of temperature and energy is represented in Fig. [9] (see an 
enlargement in Fig. [TUl) . Its description is similar to that of Sec 14. 3[ For M > Mumit, or 
equivalently n > fJ^c = '^^2 — 22.4. there exists a minimum energy Ef. = —AcGM^/R"^ and a 
minimum temperature Tc = GM/ {rjcR^) below which there is no equilibrium state. 

15 I ^ ^ ^ 
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Figure 8: The mass-energy relation for white dwarf stars (T = 0) in D = 4. Self-confined 
white dwarf stars are marginally stable. They have a unique mass Mumu independent on their 
radius. For M < Mumit, the white dwarf star is box-confined. There is no equilibrium state 
with M > Miimit- 



4.5 The case L> > 2(1 + ^2) 

The caloric curves for D > 2(1 + -\/2) are similar to those of Sees. 14.31 and 14.41 There are, 
however, two main differences. For D > 10, the classical spiral ceases to exist [14]. Thus, 
the caloric curve does not wind up as /i — > -|-oo contrary to Fig. [TOl On the other hand, for 
D > 2(1 -|- it is not possible to construct self-confined white dwarf stars [12j. This is just a 
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Figure 9: Caloric curves in D = 4 for different values of the degeneracy parameter. For 
fi > fic = 22.4, there is no equilibrium state if the temperature and the energy are too low. 
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Figure 10: Same as Fig. [9]for larger values of fi showing the developement of the classical spiral 
recovered for fi —>■ +oo. 
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mathematical curiosity since complete white dwarfs stars are unstable for D > A anyway. This 
property changes the unstable branch of the A — t]p diagram without consequence on the caloric 
curves. The A — t]p diagram is represented Figs. [11] and [121 For D > 2(1 + V2), it displays an 
infinity of cusps towards the singular solution {Agj-qp^s), see Fig. M For D = 2(1 + ^2), there 
is just one cusp (see Fig. [TT]) and the Lane-Emden equation (15^ can be solved analytically. 
This corresponds to the Z)-dimensional Schuster solution obtained for n = [12]. In that case, 
we find explicitly 

»5 = (73) 



[1 



4(2+^2) J 

The normalized mass and the normalized energy can be expressed as 

VP 



As = -2(1 + V2) 
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Figure 11: The mass-energy relation for white dwarf stars (T = 0) in D = 2(1 + V2). 



4.6 The case D = 2 

Let us now consider smaller dimensions of space. The dimension D = 2 is critical concerning 
gravitational collapse as discussed in |14j. For D = 2, the relevant Fermi integrals are /q and 
Ii. By definition, 

m = / (76) 



1 + te^ 

Changing variables to y = e^, we easily find that 



/o(t) = ln(l + i). (77) 
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Figure 12: The mass-energy relation for white dwarf stars (T = 0) in D = 5.1 



Therefore, the Fermi-Poisson equation fl20l) becomes 



^(0) =^'(0) = 0. 
On the other hand, using the identity f|T9|) . giving 



t \ t 



(78) 
(79) 

(80) 



one finds that 



r ln(l -x) , , . / In 
h{t) = - / ^ >-dx = -Li2(--), 



where Li2 is the dilogarithm. 

Consider first the Fermi gas at T = 0. In D = 2, a white dwarf star is equivalent to a 
polytrope with index 713/2 = 1. The Lane-Emden equation can then be solved analytically and 
we obtain 6 = Jo{^), where Jq is the Bessel function of zeroth order. The density drops to zero 
at ^1 = ao,i — 2.40, the first zero of Jq. Considering the mass-radius relation fISBl) in D = 2, 
we see that the radius is independant on mass. Therefore, complete white dwarf stars in two 
dimensions all have the same radius. It can be written in terms of fundamental constants as 



R* — J— 



2n \gm^G 



1/2 



0.27 



h 



3/2(71/2 • 



f82) 



The relation between the mass and the central density of the white dwarf star is 

/l2 



M 



Po 



Atc^ giv?G 



^3) 
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where d[ = Jg(ao,i) — —0.52. Thus, the density profile of a two-dimensional white dwarf star 
can be written 

pir)=PoM^-f). (84) 

This is similar to the vorticity profile of a minimum enstrophy vortex in 2D hydrodynamics 
|30[ |3T] . The energy of a complete polytrope of index n in D = 2 is E = —{n — l)GM^/8 + 
{l/2)GM'^ln{R^/R) with the convention $(i?) = [12]. Therefore the energy of a 2D white 
dwarf star is 

E = lGMHn{^). (85) 

Two-dimensional white dwarf stars exist for any mass M and they are stable. Noting that 
R^/R = (/i=K//i)^/^ = ill y/Jl where /z = Ati'^tiqGR^^ we can write the normalized energy of the 
self-confined white dwarf star as 

A = llnf^Y (86) 



2 V6 / 

Let us now consider the case of incomplete white dwarf stars that are confined by the box 
{R^ > R). This corresponds to /i < ^f. Using Eq. fl52|) . we find that a = ^/JI. Then, using the 
results of [12], we find that the normalized energy of a box-confined white dwarf star in two 
dimensions is 

A = -i (87) 

We now consider the self-gravitating Fermi gas at finite temperature T 7^ 0. According to 
Eq. ( |22l) we have a = ^/Jl. Using Eq. (1241) . we obtain 

T] = j3GMm = ^ ^'{y/JI). (88) 

We need to calculate the energy specifically because the expression ( |32|) breaks down in D = 2. 
The kinetic energy K = J pd^r can be written 

K 1 

h{ke^)idi. (89) 



GM2 7^2 

On the other hand, using an integration by parts, the potential energy is given by 

W = / {V^fdh, (90) 



where we have taken $(/?) = 0. Introducing the dimensionless quantities defined in Sec. 12. 2[ 

we get 

W 1 

^'{ifidi. (91) 



GM2 2^72 

Summing Eqs. ( 189|) and (l9Ti) . the total normalized energy of the Fermi gas in two dimensions 
is 

"^-wr^^ '''' 

The corresponding caloric curve is plotted in Fig. [131 For /i +00, we recover the classical 
caloric curve displaying a critical temperature ksTc = GMm/4 [l3j. Below T^, a classical gas 
experiences a gravitational collapse and develops a Dirac peak [H]. When quantum mechanics 
is taken into account, the collapse stops when the system becomes degenerate. The Dirac 
peak is replaced by a fermion ball surrounded by a dilute halo. At T = 0, we have a pure 
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Fermi condensate without halo. This is the ground state of the self-gravitating Fermi gas 
corresponding to the vertical asymptotes in Fig. [131 For H < Q (incomplete white dwarf stars), 
the minimum energy is given by Eq. (IHTjl and for /i < (complete white dwarf stars) by Eq. 
( 186]) . This discussion concerning the difference between Dirac peaks (for classical particles) and 
fermion balls (for quantum particles) in the canonical ensemble remains valid for 2 < D < 4. 
Note also that there is no collapse (gravothermal catastrophe) in the microcanonical ensemble 
in D = 2 [32l[Tl]. 
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Figure 13: Caloric curve in D = 2 for different values of the degeneracy parameter: = 
2, 10, 100, 10^, 10^, 10^, 10®. For /i +oo, we recover the classical caloric curve displaying a 
critical temperature Tc. Below Tc the system is expected to collapse and create a Dirac peak 
("black hole"). When quantum mechanics is accounted for, the "black hole" is replaced by a 
"fermion ball" . This result is generally valid for 2 < D < A. 



4.7 The case D <2 

We finally conclude by the case D < 2 (specifically D = 1). First, we consider the Fermi gas 
at T = 0. The A — t]p curve which gives the energy of the star as a function of its mass is 
represented in Fig. [TH Since the curve does not present turning points, all the white dwarf star 
configurations are stable. According to Eq. for _D < 2, the mass M of a complete white 

dwarf star increases with its radius R^:. Therefore, for M < M^{R) the star is self-confined and 
for M > M^:{R) it is restricted by the box. There exists a stable equilibrium state at T = 
for all mass. In terms of the dimensionless mass rjp, complete 71,3/2 polytropes correspond to 
Vp ^ and incomplete 71,3/2 polytropes to rjp > uJn/2- This situation is reversed with respect 
to that of Fig. d 

The caloric curve for arbitrary temperature and energy is represented in Fig. [T3 For 
H +00 (i.e. ^ — >■ 0), we recover the curve obtained in [13] for classical isothermal systems. 
The caloric curve f3{E) is monotonic. Therefore, there is no phase transition for D < 2. Thus, 
the change in the caloric curve due to quantum mechanics is not very important since an 
equilibrium state (global maximum of entropy or free energy) already exists for any accessible 
energy E and any temperature T in classical mechanics. Quantum mechanics, however, changes 
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Figure 14: The mass-energy relation for white dwarf stars (T = 0) in D < 2 (specifically 
D = 1). There exists an equilibrium state for all mass. The white dwarf star is self-confined if 
M < M^{R) and box-confined if M > M^{R). 




Figure 15: Caloric curve in D = 1 for different values of the degeneracy parameter (various 
system sizes). 
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the ground state of the system. The ground state of the self-gravitating Fermi gas (T = 0) 
corresponds to a white dwarf star configuration. Its structure (radius, energy) is determined 
by the intersection between the A — rjp curve in Fig. [THand the hne defined by Eq. fl62|) . The 
"white dwarf" is complete {R^, < R) for /i > ^^{D) and incomplete (i?* > R) otherwise. For 
/i > fi^:{D), the normalized energy of the white dwarf is given by Eq. (jSl])- This is the ground 
state of the self-gravitating Fermi gas corresponding to the asymptote in Fig. [T51 In D = 1, it 
is possible to obtain more explicit results. Using the results of |12], for 77,3/2 = 1/2 polytropes, 
we have ^1 = {37r/Ay/^T{5/3)/T{7/6) ~ 1.49 and \e[\ = 2/^3 ~ 1.15. Therefore, uy2 = 0.349 
and /i* = 0.846. For ii > = 0.846, the normalized energy of a complete white dwarf star 
(ground state) is 



5 Conclusion 

In this paper, we have studied how the dimension of space affects the nature of phase transitions 
in the self-gravitating Fermi gas. Since this model has a fundamental interest in astrophysics 
[H] and statistical mechanics [5], it is important to explore its properties thoroughly even if we 
sacrifice for practical applications. It is well-known in statistical mechanics that the dimension 
of space plays a crucial role in the problem of phase transitions. For example, concerning the 
Ising model, the behaviour m. D = 1 and D > 2 is radically different [33]. We have reached 
a similar conclusion for the self-gravitating Fermi gas. The solution of the problem 'm. D < 2 
does not yield any phase transition. In D = 2, phase transitions appear in the canonical 
ensemble but not in the microcanonical ensemble. In D > 2, phase transitions appear both in 
microcanonical and canonical ensembles in association with gravitational collapse. The beauty 
of self-gravitating systems, and other systems with long-range interactions, is their simplicity 
since the mean-field approximation is exact in any dimension. Therefore, the mean-field theory 
does not predict any phase transition for the self-gravitating Fermi gas in Z) = 1, contrary to 
the Ising model. 

At a more philosophical level, several scientists have examined the role played by the di- 
mension of space in determining the form of the laws of physics. This question goes back to 
Ptolemy who argues in his treatise On dimensionality that no more than three spatial dimen- 
sions are possible in Nature. In the 18*^ century, Kant realizes the deep connection between the 
inverse square law of gravitation and the existence of three spatial dimensions. In the twentieth 
century, Ehrenfest [13] argues that planetary orbits, atoms and molecules would be unstable in 
a space of dimension Z) > 4. Other investigations on dimensionality are reviewed in the paper 
of Barrow [31] . Although we ignored this literature at the begining, our study clearly enters in 
this type of investigations. We have found that the self-gravitating Fermi gas possesses a rich 
structure and displays several characteristic dimensions D = 2, D = A, D = 2{1 + \/2) and 
D = 10. Moreover, as already noted in [12], the dimension D = 4 is critical because at that 
dimension quantum mechanics cannot stabilize matter against gravitational collapse, contrary 
to the situation in D = 3. Interestingly, this result is similar to that of Ehrenfest although it 
applies to white dwarf stars instead of atoms. The dimension D = 2 is also critical as found in 
|14] and in different domains of physics. Therefore, the dimension of our (macroscopic) universe 
D = 3 plays a very special role regarding the laws of physics (this is illustrated in Fig. [T]). 
Following the far reaching intuition of Kant, we can wonder whether the three space dimensions 
are a consequence of Newton's inverse square law, rather than the opposite. We note also that 
extra-dimensions can appear at the microscale, an idea originating from Kaluza-Klein theory. 




(93) 
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This idea took a renaissance in modern theories of grand unification. Our approach shows that 
already at a simple level, the coupling between Newton's equations (gravitation) and Fermi- 
Dirac statistics (quantum mechanics) reveals a rich structure as a function of D. Relativistic 
effects will be considered in a forthcoming paper [29] . 

Finally, our study can shed light on the mathematical properties of the Vlasov-Poisson sys- 
tem. Indeed, there is a close connexion between collisionless stellar systems and self-gravitating 
fermions [191 [351 [20l |2T] . For example, the fact that the Vlasov equation does not blow up (i.e., 
experiences gravitational collapse) in D = 3 for non singular initial conditions can be related 
to a sort of exclusion principle, as in quantum mechanics. Due to the Liouville theorem in 
/i-space, the distribution function must remain smaller that its maximum initial value f < r]o 
and this prevents complete collapse [Sni [36] , unlike for collisional stellar systems |9] described by 
the Landau-Poisson system. Since quantum mechanics cannot arrest gravitational collapse in 
D > 4 (for sufficiently low energies), this suggests that the Vlasov-Poisson system can probably 
blow up for D > 4. This remark could be of interest for mathematicians. 



Acknowledgements I acknowledge interesting discussions with P. Biler, T. Nadzieja and 
R. Stahczy. I also thank B. Dougot for encouragements. 



25 



References 

[1] Dynamics and thermodynamics of systems with long range interactions, edited by Dauxois, T, 
Ruffo, S., Arimondo, E. and Wilkens, M. Lecture Notes in Physics, Springer (2002). 

[2] T. Padmanabhan, Phys. Rep. 188, 285 (1990). 

[3] P. Hertel and W. Thirring, in: Quanten und Felder, edited by H.P. Diirr (Vieweg, Braunschweig, 
1971). 

[4] P.H. Chavanis, Phys. Rev. E 65, 056123 (2002). 

[5] P.H. Chavanis, in Dynamics and thermodynamics of systems with long range interactions, edited 
by Dauxois, T, Ruffo, S., Arimondo, E. and Wilkens, M. Lecture Notes in Physics, Springer 
(2002) ||| cond-mat70212223| . 

[6] P.H. Chavanis, in: Proceedings of the Fourth International Heidelberg Conference on Dark Matter 
in Astro and Particle Physics, edited by Klapdor-Kleingrothaus, H.V. (Springer, New- York, 2002) 
astro- ph/0205426| . 

[7] P.H. Chavanis, Astron. Astrophys. 381, 340 (2002). 

[8] V.A. Antonov, Vest. Leningr. Cos. Univ. 7, 135 (1962). 

[9] D. Lynden-Bell and R. Wood, Mon. Not. R. Astron. Soc. 138, 495 (1968). 

[10] S. Chandrasekhar, An Introduction to the Theory of Stellar Structure (Dover, 1942) 

[11] R.H. Fowler, Mon. Not. R. Astron. Soc. 87, 114 (1926). 

[12] P.H. Chavanis and C. Sire, Phys. Rev. E 69, 016116 (2004). 

[13] P. Ehrenfest, Proc. Amst. Acad. 20, 200 (1917). 

[14] C. Sire and P.H. Chavanis, Phys. Rev. E 66, 046133 (2002). 

[15] S. Chandrasekhar, Astrophys. J. 74, 81 (1930). 

[16] P.H. Chavanis, preprint 

[17] H.J. de Vega and N. Sanchez, Nucl. Phys. B 625, 409 (2002). 

[18] R. Stahczy, preprint 

[19] D. Lynden-Bell, Mon. Not. R. Astron. Soc. 136, 101 (1967). 

[20] P.H. Chavanis and J. Sommeria, Mon. Not. R. Astron. Soc. 296, 569 (1998). 

[21] P.H. Chavanis, Statistical mechanics of violent relaxation in stellar systems. In: Proceed- 
ings of the Conference on Multiscale Problems in Science and Technology (Springer, 2002) 
[astro-ph/0212205j . 

[22] J. Katz, Mon. Not. R. Astron. Soc. 183, 765 (1978). 

[23] J. Katz, Found. Phys. 33, 223 (2003). 

[24] P.H. Chavanis and M. Rieutord, 2003, Astron. Astrophys. 412, 1 (2003). 



26 



[25] P.H. Chavanis, Astron. Astrophys. 401, 15 (2003). 

[26] J. Katz and 1. Okamoto, Mon. Not. R. astr. Soc. 317, 163 (2000). 

[27] P.H. Chavanis and I. Ispolatov, Phys. Rev. E 66, 036109 (2002). 

[28] B. Stahl, M.K.H. Kiessling and K. Schindler, Plane. Space Sci. 43, 271 (1994). 

[29] P.H. Chavanis, in preparation 

[30] P.H. Chavanis and J. Sommeria, J. Fluid Mech. 356, 259 (1998). 

[31] P.H. Chavanis, Phys. Rev. E 68, 036108 (2003). 

[32] J. Katz and D. Lynden-Bell, Mon. Not. R. Astron. Soc. 184, 709 (1978). 

[33] K. Huang, Statistical mechanics (Wiley, New- York, 1987) 

[34] J.D. Barrow, Phil. Trans. E. Soc. Lond. A 310, 337 (1983). 

[35] P.H. Chavanis, J. Sommeria and R. Robert, Astrophys. J. 471, 385 (1996). 

[36] R. Robert, Class. Quantum Grav. 15, 3827 (1998). 



27 



